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^ ! Abstract 

OS : 

. In this paper we investigate some properties, including causality, of a particular class of relativistic dissipative 

fluid theories of divergence type. This set is defined as those theories coming from a statistical description of matter, 
in the sense that the three tensor fields appearing in the theory can be expressed as the three first momenta of a 
suitable distribution function. In this set of theories the causality condition for the resulting system of hyperbolic 
partial differential equations is very simple and allow to identify a subclass of manifestly causal theories, which are 
so for all states outside equilibrium for which the theory preserves this statistical interpretation condition. This 
subclass includes the usual equilibrium distributions, namely Boltzmann, Bose or Fermi distributions, according 
to the statistics used, suitably generalized outside equilibrium. Therefore this gives a simple proof that they are 
causal in a neighborhood of equilibrium. We also find a bigger set of dissipative divergence type theories which 
. are only pseudo-statistical, in the sense that the third rank tensor of the fluid theory has the symmetry and trace 

properties of a third momentum of an statistical distribution, but the energy-momentum tensor, while having the 
^*) , form of a second momentum distribution, it is so for a different distribution function. This set also contains a 

■ subclass (including the one already mentioned) of manifestly causal theories. 

O 

!£) 1 Introduction 

■ 

O i In the last few years there has been a considerable effort to understand dissipation in relativistic theories of fluids. 

Straight forward generalizations to relativity of the Navier-Stokes scheme Q resulted in systems with an ill-posed 
!h i initial value formulation^]. Even if these generalizations would have worked, they would have yielded a parabolic 
^P' system, well-posed in the mathematical sense, but unacceptable on physical grounds due to the presense of infinite 
>• propagation velocities. Thus, alternative theories were proposed resulting in a formalism having an extended number 
of dynamical variables and where dissipation is at the microscopic level completely different to the standard parabolic 
dissipation of Navier-Stokes equations, but which behaves at measurable scales in the same way Q, Q. Due to the 
C3 problems encountered on the generalizations alluded above, one of the basic requirements to be checked on these 
alternative theories was their relativistic causality || , || . 

Among the alternative theories, we specialize to the case of theories of divergence type ||, j7j, that is theories with 
the following structure 

V Q A^ a = (1) 
V a T ab = (2) 

VaA abc = jbc (3) 

where the tensors A abc (tensor of fluxes) and I bc (dissipation-source tensor) are supposed to be functions on a smaller 
set of fields, which in principle can be taken to be the direct observables of the theory, namely the particle number N a 
and the energy-momentum T ah of the fluid. The tensor I ab represents the non-equilibrium interactions in the fluid, 
and vanishes at the local equilibrium fluid configurations, i.e. at the momentarily static configurations. We focus on 
this particular type of theories, not just because their structure is simpler, but also because: a) They already contain 
enought degrees of freedom to account for the usual description of nonrelativistic dissipative fluids, b) Since these 
theories generically develop shock waves, it is desiderable their structure has this form in order for the shock wave 
solutions to make sense. 
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The requirement of an entropy law, that is the existence of an extra vector field whose divergence, by the sole 
virtue of the above equations, is a function of the basic fields at the point (and not of any of its derivatives), puts 
severe restrictions in the theory. In particular it implies the existence of a single generating function x, and variables 
(C,C a ;C afc )) with ( negative, C a future directed timelike, and ( ab symmetric and trace free, such that 

N > = ik (4) 



d 2 x 

d( a d( b 

d 2 x 
d( a d( b 



(6) 



That is, this generating function determines completely the principal part, in the sense of the theory of partial 
differential equations, of the above system, and so its causal properties. 

Non-dissipative fluids can also be described with divergence type theories by simply choosing a generating function 
of the form x(CC a )i an d setting to zero the tensor I ab . This kind of generating function describes a perfect fluid 
theory, with x related with the fluid's equation of state, as explained in section II and Appendix A. 

In this paper we first define a subclass of non-dissipative divergence type theories that we call statistical. This 
perfect fluid theories can be characterized as those having a statistical origin, in the sense that the particle-number 
current and the stress-energy tensor can be expressed as suitable linear combination of the first two momenta of some 
distribution function. This requirement singles out a generating function for the perfect fluid theory as a functional 
of the corresponding distribution function, which determines completely the dynamics, thus in particular its causal 
properties. We obtain a very simple sufficient condition on the associated distribution function, to ensure causality of 
the whole perfect fluid theory. We call the theories satisfying these condition manifestly causal theories. They include 
as particular case the Boltzmann, Bose, and Fermi gasses. 

From the statistical non-dissipative divergence theory we can easily define a dissipative divergence theory, keeping 
its statistical origin. We shall say in this case that the dissipative theory thus obtained is a natural extension of the non- 
dissipative one. This extension, has an important property: If the original non-dissipative theory is manifestly causal, 
then its dissipative extension is also manifestly causal. In particular, since the theories of equilibrium Boltzmann, 
Bose, and Fermi gasses can be cast in the form of manifestly causal non-dissipative theories, we conclude their natural 
extensions are also manifestly causal, thus considerably generalizing and simplifying of the works of [|j , |J , JlO| , JTTJ . 
In spite of its naturality, this extension has some drawbacks, as explained in Section II, it is only well defined for some 
values of the dissipation variables, which for instance do not cover any whole neighborhood of the local equilibrium 
submanifold. Thus, some compatible, but otherwise arbitrary, extension has to be made in order to cover the remaining 
values the dissipation variables can take. 

Later on, in section III, we shall relax the statistical condition and allow theories which have a pseudo- statistical 
origin, that is where the three tensors fields, N a , T ab , and A abc are averages over a future mass shell, but for different 
(although related) distribution functions. Again we shall find for this case a subclass of manifestly causal theories 
which includes the former one, and again show that manifestly causal non-dissipative theories have natural dissipative 
extensions which are also manifestly causal. 

In section IV we shall discuss the expression for the entropies of this theories. The requirement of an entropy law is 
so strong that even singles out a unique form of the entropy as a functional of the generating function. The entropies 
we deal here are dynamical entropies in the sense that they appear as further relations (equations) between the fields 
as a consequence of its dynamics. We find that, in contrast with the usual definitions of equilibrium entropies, for 
the statistical theories they are linear functionals of the distribution functions. Nevertheless, when evaluated at given 
equilibrium distribution functions, e.i. Boltzmann, Bose or Fermi distribution functions, they coincide with the usual 
formulae coming from sum-of-states considerations. 

Finally, for completeness, we include two appendices, one is a short review of the relativistic dissipative theories 
of divergence form, and the other of relativistic statistical mechanics. 



2 Statistical Theories 

The divergence type fluid theories are briefly reviewed in Appendix A. The simplest particular case of these theories is 
when the generating function x does not depend on dissipative variables, and therefore is non-dissipative. In this case 
the theory describes a perfect fluid, with the particle-number current and the stress-energy tensor having the usual 
structure, namely, 



N a 



A* 
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where q = g + u a u , and £ a = p,u a , with u a u a = — 1. So we have the particle-number of the fluid n — —x,(v> the 
energy density p = X,/j.h an( ^ pressure p = — ^f- in terms of x- It can be seen that in this theory p is one over the 
temperature of the fluid, and £ is a chemical potential per unit temperature. The choice x = x(C>C a ) implies that 
the tensor A abc is identically zero, and this is consistent with the choice I ab = 0, for the dissipation-source tensor. In 
these theories the expression for the entropy current given in the Appendix reduces to the usual one, S a = nsu a , with 
s = (p+p)/(nT) - C, and it satisfies V a 5 Q = 0. 

We say that a perfect fluid theory is of statistical type, if its generating function x(C;C a ) can be- written in the 
following way 

X((,C)= [ f(( + PaC)duj (7) 



where the integral is on the future mass shell p a p a — —m 2 , and / is a smooth scalar function on the negative real line, 
which for large (negative) values decays fast enought as to make the integral well defined and x(C,( a ) smooth in all 
its variables. Notice that £ < 0, since it is a chemical potential, while p a C a < 0, since it is the scalar product of two 
future directed time-like vectors. 

The simplest example is when f(x) = k 2 e x / k , where k is Boltzmann's constant, namely Boltzmann's gas. Indeed, 
the expressions for the particle-number current and the stress-energy tensors are in this case given by 

N a = J Pa e «+P»C)/fc dulj 

Tab = I PaPb e«+^ a )/ fc dw. 



That is, we interpret N a and T a b as the two first momenta of the Boltzmann's equilibrium distribution function. 

For an arbitrary statistical perfect fluid it is easy to see that the x defined as in (Jfy corresponds to a statistical 
theory with a distribution function given by /". 

An interesting property of this integral representation is that causality is easy to analyze. Indeed following || (see 
Appendix A) , we say that a perfect fluid theory is causal if: 

t a E a = l -t a M a AB Z A Z B < 

for all perturbations Z A = (<5£, 5( a ) and all future-directed timelike vectors t a . But for statistical perfect fluid theories, 
using the expression of x m terms of mass shell integrals, this condition becomes: 

taE a = \l {t a p a )f" (S( + PadCf , (8) 

for all perturbations (5(,8( a ). Thus we can isolate from the set of all statistical theories a subclass of manifestly 
causal theories, namely those having /"' ^ 0, /"' > 0. [] Notice that Boltzmann's gas is inside this subclass, which 
also includes Bose and Fermi's gases, since for them we have 

(according of the value of e, 0, 1 or —1, we have respectively Boltzmann, Bose or Fermi gases), and it is easy to see 
that for all of them /"' > 0. The generating functions x f° r Bose and Fermi distribution functions are obtained by 
straightforward integration. 

We now turn to the problem of extending these theories outside equilibrium. 

The extension outside equilibrium we propose is the following: Given a statistical perfect fluid theory characterized 
by a function /, we define a dissipative divergence theory of statistical type by the generating function 

X(C, C, C b ) = J f(C + Pa? + PaPbC") dw. (9) 

Note that this extension is only valid for values of C ab such that the argument of / is negative, while in general 
( ab can take values for which the argument is positive q. Since £ and p a ? are already negative for all p a in the future 
mass shell, the values of £ ab for which the argument is negative for all values of the momentum in the future mass 
shell form a cone, C~ = {C afc | C ab = C fca , 9ab? h = 0, l a lb? b < 0, V l a null}, of maximal dimension. Thus this 
extension is not unique, the equilibrium behavior of the statistical fluid only gives information encoded in an / which 
is only defined for negative values of the real line, while here we seem to need an / defined on the whole line. But the 
situation here is even worse, as shown in the following lemma. 

1 Note that there are theories which are not manifestly causal, but nevertheless are causal. 

2 Near equilibrium the theory should behave, at least in some respects, like Eckart's theory. But in this theory there is a relation between 
£ ab and the derivatives of the flow velocity which renders totally unphysical any imposition on the sign of £ ab . 
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Lemma 1 The function, 

F{c) = I f{-x + cx 2 )x 2 dx, 



where f is any smooth, positive definite function, which is of compact support or decays faster than x 4 is discontinuous 
at c = 0, having there a finite limiting value from the left and diverging from the right. 

Proof 

To prove that the limiting value from the left is finite just note that if c is negative, then the argument of / just 
grows in absolute value and so the decay assumptions on / imply convergence. To prove the statement about the the 
limiting value from the left, take a double step function (a small square), sc(x), smaller than /, then we have, 

/>oo 

F(c) > G{c) = / sc(-x + cx 2 )x 2 dx. 



So it suffices to prove the divergence on the limit for the function G(c). Assume for simplicity that sc(x) is different 
from zero only in the interval [—1, 1], and that there its value is 1. For positive values of c, the integral above has 
contributions from two intervals. If c is small, those intervals can be calculated up to first order in c obtaining [1 + c, 0] 
and [- + 1, - — 1]. The contribution to the integral from the first one is finite for all values of c, while the second can 
be estimated using the mean value theorem to be bigger than ^ + 1 + 0(c) and so we see that it diverges for c going 
to zero from positive values. 

Q.E.D. 

This Lemma gives a strong argument against the possibility of extending smoothly the definition of x to values 
of ( ab such that p a PbC ab > 0, by simply extending the definition of / to positive values, no matter how smooth or 
how strong a decay condition we impose on the extension. But of course, there exist a lot of smooth extensions of the 
generating function x to the presently forbidden values of ( ab . Any of those, essentially arbitrary, smooth extensions 
will be assumed to have been made in what follows, in particular we require that extension in a neighborhood of 
equilibrium, that is, in a neighborhood of the apex of the cone C'~ . The results on causality near equilibrium do not 
depend on the particular extension chosen outside . 

For these statistically extended dissipative fluids theories also a simple expression for the causality condition is 
easy to obtain, namely: 

t a E a = \J {t aP a )f" (5( + pJC + PaPbSC") 2 
So we have have the following result: 

Theorem 1 Let function f : R~ — > R defining a statistical perfect fluid be C 3 , and such that the equilibrium theory 
is well defined. If the equilibrium theory is manifestly causal, i.e. f" > 0, then the extended dissipative theory is also 
causal in a neighborhood of equilibrium. 

Proof 

The above expression shows causality for all values of ((, ( a , ( ab ) such that ( ab 6 C~ . Since the cone is of maximal 
dimension, partial derivatives along directions inside the cone suffices to determine completely the differentials of x a t 
the apex of the cone, that is at equilibrium. Thus the smooth extension outside the cone can not change the causality 
properties of the equilibrium configuration. The result extends to a neighborhood of equilibrium trivially by noticing 
that in our setting causality is a continuous property. 

Q.E.D. 

Particular examples of this theorem are the dissipative extensions of Boltzmann, Bose and Fermi's gases, given by 

1 



f" 



e -«+PaC a +PaPb< ab )/k 



according of the value of e is 0, 1 or —1 respectively. In all cases is easy to see that /"' > 0, so they are smooth and 
manifestly causal for all values of the parameters for which the integral expression converges, in and off equilibrium. 
This generalize and simplify our previous work |p"lf , which in turn was a generalization of other results J§| , jl0| . 

3 Pseudostatistical Theories 

There is an even larger set of theories defined in terms of the integral of certain function / on a future mass shell, for 
which it is straight forward to find sufficient conditions on / such that the resulting divergence theory is causal. We 
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say that a dissipative fluid of divergence form is of pseudo-statistical origin, if its generating function x can be written 
as 

^CCCH / f(C+PaPbC b ,PaC)dLU 

where the integral is on the future mass shell, and / = /(&, y) is a scalar function of two variables, x = C + PaPbC ab 
and y = p a C, a . We call them pseudo-statistical because of the tensors N a and A ahc can be thought as coming from a 
distribution function f xy , ^] while the stress-energy tensor T ab can be thought as coming from a different distribution 
function, f yy , in fact, 



N a = J p a f xy duo 

Tab = / PaPb fyy duJ 



Aabc = J PaP<bP C fxy did . 

Here the values of the variables (CO have the same restrictions as in section above. The lemma 1 does not apply 
here and there are cases, like Boltzmann, where one can obtain a pseudo-statistical extension for all values of C, ab 
starting from an equilibrium statistical theory. 

The pseudo-statistical theories are interesting since again the causality condition is simple and it is easy to impose 
a sufficient condition such that the theory is manifestly causal. Indeed, it is easy to see that in this case the causality 
condition is, 

t a E a = \ J(t aP a ) [fxxy(SC) 2 + fyyy (PbSf) 2 + fxxy (p<dPe>S( de ) 2 + 

2/xvA (PbSC b ) + VxxyK (p<dPe>5( de ) + 2f xyy ( Pb 8C, b ) {p <d p e> 5C, de )} diO . 

Rearranging terms we have 

t aE a = \ J (taP a ) [fxxy(S(+P <d Pe>SC de ) 2 +f vyy (pbSC b ) 2 +2f xyy (SC+ P<dPe> SC: de ) (rf) 

and defining Sx = S( + p< a Pb>o~( ab and Sy — p a o~( a , 



duj, 



taE a = \ I {t a p a ) [fxx y (5x) 2 +f yyy (5y) 2 + 2f xyy Sx5y] dw. (10) 



Thus, if we assume f xxy > 0, f yyy > and (f X y V ) 2 < f y yyf xxy , we obtain, 



\ ( (taP") (yi^V \5X\ ~ yfjyy'y' \5y\j (fa), 



t a E a < 
~ 2 

and again isolate a subclass of dissipative theories (containing the one previously described) which are manifestly 
causal, namely those which satisfy, 

fxxy ^ fyyy — ^ (.fxyy^) — fyyyfxxy 

We thus have an obvious generalization of the previous theorem. 

Theorem 2 Let function f : RT x Rr — > R defining a perfect (equilibrium) pseudo statistical fluid be C 3 , and such 
that the equilibrium theory is well defined. If the equilibrium theory is manifestly causal, then the extended dissipative 
theory is also causal in a neighborhood of equilibrium. 

Remark: There are causal /'s whose derivatives does not satisfies conditions above. Even more, if we extend the 
theory in a pseudostatistical way for positive arguments, then the decay properties would in general imply the above 
condition can not be imposed, so even in this case, causality for naturally extended theories seems to hold only in a 
restricted domain of possible field variables. 

The generating functions x associated with a statistical or a pseudo-statistical dissipative divergence type theories, 
can be thought as particular solution of equation ( |30|) of Appendix B. This equation is obtained by imposing some 
conditions on the divergence type theories, related with the symmetry and trace of the third rank tensor, A abc , of the 
theories. Those conditions come naturally from kinetic theory and are explained in that Appendix. From the above 
subclass of theories it seems then that there arc globally defined, and nicely behaved solutions to those equations, but 
they are not statistical, as one would have liked them to be, they are at most pseudostatistical. 

3 Here the subindices indicate partial derivatives with respect to the argument signaled by the index. 



■5 



4 Entropy 



As seen in the introduction and Appendix A, the requirement of an entropy law is a very strong condition which 
substantially constraints the kinematics and dynamics of the fluid theories, since imposes that it can only depend on a 
freely given generating function, \, (essentially an equation of state), and a -basically freely given- dissipation source 
tensor I ab . Even more, it gives a fixed formula for the entropy as a functional of the arbitrary generating function 
defining particular theories. This entropy is of dynamical origin and in principle it is not related to the entropy concept 
coming from information theory as applied to equilibrium configurations. In fact, the entropy law property can be 
divided in two conditions the more stringent one is the existence of an extra vector field whose divergence is, as a 
consequence of the other equations, a local function of the dynamical field (and not of their derivatives) , and a weaker 
one which requires this last function to be non-negative and to vanish at equilibrium. The first condition is the one 
which restricts the theory, the one which gives its rigidity, and is the one which determines completely the entropy 
functional, while the other condition, in this picture, just puts mild conditions on the dissipation-source tensor. We 
shall see in what follows how these conditions work for the natural extensions of statistical theories. 
In the theories of divergence type the entropy is given by, (see Appendix A.) 

S a = - (N a - ( b T ab - ( bc A abc 



where ( A = ((,( a ( ab )- For the naturally extended statistical theories this translate into, (See Appendix B.) 

Sa= J Pa (/' - C A P A f") duj 

where P A = (l,p a ,p a p b ). 

We observe that this entropy is a linear functional of the distribution function, and so there is no convexity property 
a this stage. Its form does not depends on any statistical counting of states assumption, once the distribution function 
/" is given, the entropy is determined. 

The relation with the usual concepts and formulae comes about when one identifies, via Boltzmann equation, the 
tensor I ab with a collision integral. That is we pretend that the dependence of this tensor with ( A is only throught 
the distribution function. In that case, the weaker condition we had in the dynamical picture, i.e. the vanishing of 
the divergence of the entropy current to second order in the dissipation variables implies, via the collision functional, 
that the global equilibrium states have the Gibbs form, which are thus linked to the statistics considerations used to 
deduce the form of that functional. To see how this works we assume the above setting and compute the entropy for 
the global equilibrium states of the usual statistics, (they enter only in the choice of collision functional used to define 
I ab ), that is the examples considered above. Thus assume that the associated distribution function to the fluid theory 
is of the form 

f" = F = ^r—, 

e A — e 

where n is a constant proportional to 1/h 3 with h the Plank's constant and for simplicity we denote X = —PA( A /k — 
— (x + y)/k = —(( + p a ( a + p a Pb( ab )/k. It is direct to see that 

f' = -^ [ln(e*-e)-X]. 

If we define A = 1 + e F/rj, then we can write e x — (i]A)/F; and 

/' = ^ ln(A). 



So the entropy current density functional S a can be written as 

S a = J Pa (f + kX f") (Lj 



k J Pa ln(A)+ln^) Fj dw 

~k J p a 



F In (^j - 5 A ln(A)^ du>. 

The last expression correspond to the usual entropy functional defined in relativistic kinetic theory whose equilibrium 
states corresponds to Bose or Fermi's equilibrium distribution function if e is 1 or —1, respectively. The case e = 
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corresponds to the Boltzmann's distribution function, and is easy to see that S a given by (0) has the usual form 



Thus the usual formulae for the entropy functionals is recuperated, but they are only valid for global equilibrium 
distribution functions, that is the properties usually assigned to entropies are only relations among certain distribution 
functions, and do not seems to hold outside that set, even at local equilibrium states. 

5 Conclusions 

We have presented a variety of dissipative divergence type theories with a statistical origin, in the sense that the tensors 
of the theory can be expressed as appropriate functions of the three first momenta of a suitable distribution function. 
This represents a relation with kinetic theory, which is manifest in the integral expression for the resulting generating 
function \- From this integral expression we could easily derive a simple condition on the associated distribution 
function for the resulting theory to be causal, even for some states even far momentarily equilibrium states. This 
condition had been easily verified for the natural extensions to non-equilibrium states of equilibrium distribution 
functions associated with Boltzmann, Bose and Fermi's statistics. The results obtained for these particular examples 
represent a great simplification of our previous work Jllj , which was a generalization of former works || , |1C[ ] . 

The dynamical entropy defined in the divergence type theories has not been related, in principle, to the entropy 
concept coming from a statistical theory as applied to equilibrium distribution functions. This dynamical entropy is 
just a vector field (constructed from the fluid fields) whose divergence is a pointwise function in the fields and its 
definition does not use anything about equilibrium configurations. This expression has a surprising form, for in the 
theories here considered it is linear on the associated distribution function. A remarkable fact for this dynamical 
entropy is that for momentarily equilibrium configurations corresponding to the distribution functions associated with 
the usual Boltzmann, Bose and Fermi statistics; it takes the familiar form encounters in statistical mechanics when 
applied to equilibrium configurations. Thus, from this point of view, the expressions usually used are only valid for a 
very restricted set of distribution functions and when so, only for some regions of the non-equilibrium configuration 
variables manifold, namely that where the natural extension holds. 

There exists a serious limitation for this integral representation of the statistical theories, as has been shown with 
the Lemma 1. This Lemma says that we can not extend smoothly the definition of x, to values of ( ab such that 
( ab PaPb > 0, simply by extending the definition of / to such values. In other words, the statistical interpretation 
of x can n °t be extended. We can, of course, extend smoothly the definition of \ to the forbidden values of ( ab 
in an almost arbitrary way. All those extensions will be, by continuity, causal in a neighborhood of the region 
where p a Pb(, ah is non-negative, in particular in a neighborhood of equilibrium. Does the limitation in the integral 
representation, have any physical interpretation?. The condition l a lbC, ab < 0, V7° null seems to be unphysical if we 
look near equilibrium, for there the theory should resemble Eckart, and in that theory it is known that there does not 
exist any restricting condition on the dissipative variable £ ab . So the problem, if this is really a problem, should be 
in the fluid approximation, which in this case is the pretention that we can describe the evolution of a distribution 
function via Boltzmann's equation via the evolution of a finite set of variables, namely ( A = (C, CS C ab )- While this 
seems to work very well at equilibrium, it might be that it fails completely away from equilibrium and for an 
appropiate description an infinite of field variables are needed. In that case a natural extension seems to be possible, 
provided the distribution function is suitably extended for positive values of its argument. This solution is of the type, 
add more dynamical fields in order to get some sort of cutoff. Other type of possible extension is to avandon the direct 
relation between distribution momenta and dynamical variables outside equillibrium, but still retain the flavor of a 
statistical theory in the sense of the representation as integral over mass shells. This is easily implemented by changing 
the argument of the distribution function in a non-linear way, for instance: x = C + PaC a + PaPbC, ah + ^(PaPbC ab ) 2 ) ■ 
This solution is of the type, add a constant, in order to get a non-linear cutoff. Both sound interesting and perhaps 
the solution lies in between. 

It has been possible to extend the proof of causality given for the statistical theories to a bigger set of divergence 
theories. For this set, the three tensors fields, N a , T ab , and A abc are also averages over a future mass shell, but for 
different (although related) distribution functions. We do not have any application for this larger class, but believe 
it should be of relevance for describing some physical phenomena. In this case extensions seems to be possible for all 
values of the non-equilibrium variables, but such extensions in general are not causal for all such values, and so do not 
seem to be very physical. Are there special extensions which from a manifestly causal equilibrium statistical theory 
yield a causal pseudo-statistical non-equilibrium theory? 



4 It is not clear in this context what it means expressions like near equilibrium and its relation with, for instance, the the size of £ ab . 



S a = —k J p a F In I — J duj 




7 



6 Appendix A: Review of dissipative fluid theories of divergence type 

Following H , [Q , H , we define a Dissipative Fluid Theory of Divergence Type as a theory having the following three 
properties: 1) The dynamical variables can be taken to be the particle- number current iV a , and the (symmetric) 
stress-energy tensor T ab . 2) The dynamical equations are 

V Q iV Q = 0, (12) 
V a T Qb = 0, (13) 
\7 a A abc = I bc \ (14) 

where the tensors A abc (tensor of fluxes) and I bc (dissipation-source tensor) are local functions of the dynamical 
variables N a and T ab , and are trace free and symmetric in the last two indices. 3) There exists and entropy current 
S a (local function of N a an d T ab ) which, as a consequence of the dynamical equations, must satisfy 

V a S a = a (15) 

where a is some positive function of N a and T ab . 

It is proved in 0, that a theory having these three properties is determined by specifying a single scalar 
generating function x an d the tensor I ab as functions of a new set of dynamical variables £iC a :C afc (with the later 
trace free and symmetric). The dynamical equations for these variables are (|l|||), with 

= W (16) 



d 2 X 
d( a d( b 

d 2 x 
d( a d( bc 



(18) 



while the entropy current is given by 



Sa = w ~ CNa ~ c " Tab ~ (bCAabc (19) 

while its divergence is 

O = -C'lab. (20) 

It is helpful to represent the collection of dynamical variables as, C, A = (£, £ a , C ab ), an d the dissipation-source tensor 
as, I a — (0,0, lab)- Equations (|l|-||) can then be written in this notation as 

M% B V a ( B = I A (21) 

with 

M a = - ^ 
AB dCad( A dC B 

The system of equations (|2l] ) is automatically symmetric, since partial derivatives commute. We say that a symmetric 
system is hyperbolic if the vector 

E a = ~ M a AB Z A Z B 

lies to the future of some space-like three-dimensional sub-space of the tangent space, for all non- vanishing Z A . The 
system is called space-time causal if furthermore E a lies within the future light cone (i.e., if E a is a future-directed 
timelike vector), for all non- vanishing Z A . The property of hyperbolicity ensures that system (^TJ) has a well-posed 
initial-value formulation, while space-time causality ensures that no fluid signals can propagate faster than light. 
Because the underlining statistical mechanics origin of fluids, it is natural to demand causality even for its dissipative 
versions, this implies conditions on the generating function x- 

Following (| we call a state, C A j i- e - t ne value of <^ A at a given time, a momentarily equilibrium state, if at that 
time the entropy production vanishes, i.e. a = 0. In the generic case this implies C ab = 0. i.e. C„ e = (OC a ,0). An 
equilibrium state or strict equilibrium state is a time reversible state. In that case it can be seen that generically not 
only Cg = (£, C\0), and I ab (( A ) = 0, but ( is constant, and ( a is Killing. For these strictly equilibrium states the 
modulus of ( a , that is fi = \J — ( a ( a , can be associated with one over the temperature of the fluid, and its direction 
u a = C°/m with the 4- velocity of the fluid. The variable £ can be associated with a chemical potential per unit 
temperature of the fluid. 
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7 Appendix B: The relation with Kinetic Theory 

We consider a distribution of identical particles in a given, fix, spacetimc. The particles interact via short-range forces, 
idealized as point collisions. A distribution function F(x,p a ) is defined Jl2| by the statement that 

F— dT, a duj 
m 

is the number of world-lines cutting an element of 3-surface dY, a and having 4-momenta p a which terminate on a cell 
of 3-area dui on the mass shell p a p a — —m 2 . This distribution function is the solution of the relativistic transport 
equation 

p a V a F(x,p d ) =C(x,p d ) 

where the derivative is along a curve in phase space which is geodesic, starting at the point where we evaluate it and 
having as tangent vector p a . C is a collision term defined requiring that 

C{x,p d ) — \/^gd 4 x 
m 

be the number of particles in the momentum range dio around pa which are created by collisions in the 4-volume 
^—gd i x, around the point x. Given an expression for the collision term, in general a functional of F, Boltzmann's 
equation determines the dynamics of F. This collision term can not be an arbitrary function of F, it has to satisfy 
certain physical properties. Usually the following general properties are required: 

1. The form of C is consistent with 4- momentum and particle number conservation at collisions. 

2. The collision term C yields a non-negative expression for the entropy production. 

To express these two conditions in a formal way, it is necessary to relate the information on the distribution function 
with macroscopic quantities, such as particle-number current or the stress-energy tensor. To do this it is useful to 
introduce, given a dsitribution function, F the set of moments associated to this distribution, which is the following 
hierarchy of totally symmetric tensors 

J ai - a „ = J p ai - ■ -Pa n F duj, n = 0, 00 

where the integral is on the future mass shell, p a p a = —m 2 . 

These moments are not independent quantities because they satisfy the following relations 

■I; -a - ai -a, '/" '•' "'"-A, i,j = l---U i^j. (22) 

By virtue of transport equation, the (n + 1) momentum satisfies, 
with the source tensor I ai ... an defined as, 

Ioi-o„ = J Pa t - ■ -Pa n C duj. 

If we identify the first momentum with the particle- number current, that is N a = J a ; and the second momentum 
with the stress-energy tensor of the fluid, that is T a \> — J a j,; then condition 1) on the collision term can be written as 



C duj = / p a C duo = (24) 

The statistical entropy current density is defined as a functional of F as follows,0 

SJx) = -— [ <f)(F) Pa du with <f){F) = (Fln{h 3 F) - 4^ Aln ( A )" 

m J \ eh i i 

and 

h 3 

A(x,p a ) = 1 + e — F(x,p a ) 

UJ 

5 This statistical entropy is not the dynamical entropy we are using. They only coincide for equilibrium states. 
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where h the Planck's constant, u> is the spin- weight (number of available states per quantum phase-cell) and e is 1 for 
Bosons and —1 for Fermions. Then condition 2) on the collision term can be written as, 

V a S a = J c/>'(F)C{F) du >0 where 0(F) = In (j^j ■ (25) 

Finally, local equilibrium distributions (or states), which are defined requiring that for them the entropy production 
vanishes, namely V a 5 a = 0. By first property imposed on the collision term, a sufficient condition is that for this F's 

This conclude a brief review of the basics Kinetic Theory we need. 

To relate Kinetic theory with the divergence type fluid theories it is necesary to make some identifications between 
the basic thensor fields of both theories. As we did above, is direct to associate the two first moments J a and J ab with 
the particle- number current N a and the stress-energy tensor T ab . These relations do not imply any restriction on the 
generating function \, at least from the symmetry and number of linearly independent elements of these tensors. Such 
a restriction on \ comes if one whishes to relate the third momentum of a distribution, J abc , with the tensor A abc of 
the divergence type theories. This is sugested because, in analogy with the divergence type dissipative fluids the first 
three momenta satisfy 

V a J a = (26) 
V Q J afc = (27) 

Va ja(bc) = j(6c> ( 2g ) 

where the symbol ( ) means symmetrization and trace free. But this analogy is only superficial and involves a very 
stringent assumption. Indeed, the above momenta equations are not a closed system of equations, for l ab is not 
a function of some finite number of variables -for instance the previous two momenta-, but it is a function of the 
distribution function, which in general is a functional of all momenta. Thus systems are analogous only for those 
situations which can be described with a distribution function which is to a certain degree of aproximation a function 
of only the first two momenta. At present there is not convincing argument that this is ever the case, nevertheless we 
proceed here with the analogy. Since J ao b — —m 2 J a = —m 2 N a , while, A ao b = 0, and assuming that I a t, = l ao , the 
only possible relation between these tensors such that both dynamical system are equivalent is: 

2 

771 

Jabc = Aabc 4~ Na 9bc (29) 

Since J a j, c is totally symmetric, this condition imposes conditions in A aoc and iV a (such that the right hand side of ( p9| ) 
be totally symmetric), which in turn impose conditions on the generating function x- These conditions are equivalent 
to the following system of equations for \ 

d 2 \ ni 2 d 2 x 

Thus, for this identification to make sense the generating function must be a solution of (|30|). If we obtain such a 
solution, we can then describe a fluid by the corresponding divergence type theory, which thus has a statistical origin 
as given by the above identification, even away from equilibrium. A family of solutions of this system of equations is 
given by 

x(C,C a ,C ab )= J f{x,y)du, (31) 

where x = £ + p a PbC ab \ V = PaC a i the integral is on the all future-directed momenta p a , and / is assumed to behave 
asymtotically in such a way that the integral converges. With subindices in / indicating differentiation, we have 



d 2 X 



Pa P<bPc> fxy du 



d( a d( bc ~ ' FaF<blc>Jx y 

,2 



771 

Pa ( PbPc + —9bc ) fxy duj 



J PaPbPc fxy dhJ + ^~gb, 



V 

I ;)2 



2 

where p<bPc> = PbPc + ^~9bc means symmetrization and trace free. So have 

d 2 x m 2 d 2 x 



- 9b]c = / dlO P[ a p b ]Pc fxy = 
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As we have seen in the section II, there seems to be that there are not pure statistical solutions, i.e. solutions of a 
the type: f(x,y) = f(x + y), with the appropiate differentiablility conditions on \ they generate. Thus the analogy 
discussed above seems to be unjustified. 
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